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ABSTRACT 

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 (z — 0.189), in 
deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. 
Our program is designed to study the dark matter distribution in 25 massive clusters, and balances 
depth with a wide wavelength coverage, 2000-16000A, to better identify lensed systems and generate 
precise photometric redshifts. This photometric information together with the predictive strength of 
our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a 
total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile gradient, 
d log E jd log r ~ —0.6 ±0.1 (r < 160kpc). We find consistency with the standard distance-redshift 
relation for the full range spanned by the lensed images, 1.01 < z < 6.03, with the higher redshift 
sources deflected through larger angles as expected. The inner mass profile derived here is consistent 
with the results of our independent wcak-lensing analysis of wide-field Subaru images, with good 
agreement in the region of overlap (~ 0.7—1 arcmin). Combining weak and strong lensing, the 
overall mass profile is well fitted by an NFW profile with M mr = (5.37±qH ± 0.26) x l0 14 M Q /h and 

a relatively high concentration, c V i r — 8.77~^o'42 ^ 0.23, which lies above the standard c-M relation 
similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of 
other lensing clusters, te — 16 ± 2" (for z s = 2.55), so the relatively large number of lensed images 
uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH 
survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are 
X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free 
from lensing bias. 

Subject headings: dark matter, galaxies: clusters: individuals: Abell 383, galaxies: clusters: general, 
galaxies: high-redshift, gravitational lensing 
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Clusters of galaxies play a direct and fundamental role 
in testing cosmological models and in constraining the 
properties of dark matter (DM), providing unique and 
independent tests of any vi able cosmology and struc 
ture formation scenario (e.g., |Lahav et al.||'l991| |Evrard 



et al. 1120021 |Broadhurst et al.||2UU5a| ILemze et al.||2Qi 
Julio et al.|2010 ) . Their extreme virial masses mean that 



unlike individual galaxies, gas cooling is not capable of 
compressing the dark matter halo, so that cluster mass 
profiles reflect directly the thermal evolution of the DM 
and the g rowth of the cosmological density field ( |Pee-| 
bles||1985| jDuffy et al.||2010j ). The capability of clusters 
to critically examine the standard cosmological model 
is now welcomed more than ever given the unattractive 
hybrid nature of the standard ACDM model derived by 
other means. 

Simulated CDM dominated halos consistently predict 
mass profiles that steepen with radius, providing a dis- 
tinctive, fundamenta l prediction for this form of DM 



(Navarro et al. ( 1996 ); NFW). Furthermore, the degree of 
mass concentration should decline with increasing clus- 
ter mass because clusters that are more massive, collapse 
later, when the cosmological background density is lower 
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e.g.. |Bullock et al.||2001[ |Zhao et al.||2003| |Neto et al. 



2007 i Cluster lensing provides a model independent 
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means of testing these fundamental predictions. Given 
an unbiased sample of relaxed clusters with high spa- 
tial resolution, one can rigorously test these basic pre- 
dictions of the standard ACDM model and contending 
scenarios. To date, only limited progress has been made 
toward these aims given the considerable observational 
challenges of obtaining data of sufficient quality for ac- 
curate weak and strong lensing work. 

Full mass profiles spanning the weak and strong lens- 
ing regimes have been constructed for only a handful 
of clusters, involving deep HST data to reliably iden- 
tify large samples of multiple images, and high quality 
wide- field imaging for careful weak-lensing (WL) work 
(e.g.,|Gavazzi et al.||2003| |Broadhurst et al.|2"005a[ |2008[ 
Umetsu fc Broadhurst||2008| |Merten et al'||2U(7^TpDTl1 
Newman et al.|2009HCoe et al.|20101|Umetsu et al.|2010[ 



2011b||7itrin et al.|2010 ). It has become clear that the in- 
ner mass profile can be accurately obtained using several 
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2009b 


2010 


2011c 



In the case of WL the 



mass profile (jKaiser fc Squires 1993 1, but much published 
work has suffered trom a significant dilution of the lens- 
ing signal by foreground objects and cluster members, 
leading to shallow profiles with underestimated Einstein 
radii. The ability of multi-color photometry to isolate 
foreground and background with refer ence to the radial 
WL si gnal has been demonstrated by |Medezinski et ah] 
(2010), so that the WL signal is found to be higher than 
earlier work, particularly so towards the center of the 
cluster. 

The initial results from combining deep strong-lensing 
(SL) work with minimally-diluted WL analyses has led 
to intriguing results, in the sense that although the mass 
profiles are well fitted by NFW-like profiles, showing the 
continuously steepening logarithmic gradient consistent 
with the expected form for CDM dominated halos, the 
concentration of matter in these halos seems to lie above 
the mass-concentra tion relation predicted by the stan- 
dard ACDM model dGavazzi et al. 2003 ; Broadhurst et al 



2005a| |Zitrin et al.||2010| | Umetsu et al.puilbp . Lensing 



bias is an issue here tor clusters which are primarily se 
lected by their lensing properties, where the major axis of 
a cluster may be aligned preferentially close to the line 
of sig ht, boosting the projected mass density observed 
(e.g., |Hennawi et al.||2007| |Corless fc King|[2009l |Oguri 
fc Bland W||2^()9l | Sereno et al.JIMtil |Morandi et al. 



201 ip . This will usually result also in higher measured 



concentrations and larger Einstein radii (e.g., Sadeh & 



Rephaeli 2008 Meneghetti et al. 2010a), though even 



with these effects taken into account there seems to be 



some discrepancy from ACDM prediction s (Og uri et al. 
20091 |Meneghetti et al.]|2011| |Zitrin et al.||2Ulla| ). While 



existing data may not support a strong conclusion that 
the observations are in significant tension with the stan- 
dard ACDM model, it is clear that a larger X-ray selected 
sample, with minimal lensing bias and excellent SL and 
WL data, is required to evaluate the significance of these 
trends. 

Several examples of high-redshift virialized clusters 
with diffuse X-ray emission are known, where the 
highest-redshift cluster selected by X-ray means is now 
established at 



2011). 



z = 2.07 (CL J1449+0856; |Gobat et al 
The most massive of these clusters is XMMU 



J2235. 3-2557 at z = 1.39 ( TRosati et aL]|2009 ) with an 
estimated total mass of M to4 (< IMpc) = (5.9 ±1.3) x 
10 14 M Q . The existence of these clusters, as well as the 
existence of evolved galaxies at high redshift, are claimed 
to be unlikely given the predicted abundance of extreme 
perturbations of clust er sized masses in the standard 



Daddi et al. 2007 2009 Collins 



ACDM scenario (e.g. 

et al.|2009|[Jee et al.|2dl)y||Richard et al.|20Tl| , pointing' 



towards a more extended early history of growth, or a 
non-Gaussian distribution of massive perturbations. 

To shed new light on these mysteries we have em- 
barked on a major project involving galaxy clusters, 
the Cluster Lensing And Super nova survey with Hubbl e 
(CLASH). For more details see Postman et al. (2011). 
The CLASH program has been awarded 524 orbits of 
HST time to conduct a multi-cycle program that will 
couple the gravitational-lensing power of 25 massive in- 
termediate redshift galaxy clusters with HST's newly en- 
hanced panchromatic imaging capabilities (WFC3 and 
the restored ACS), in order to test structure formation 
models with unprecedented precision. The CLASH ob- 
servations, combined with our wide-field optical and X- 
ray imaging, represent a substantial advance in the qual- 
ity and quantity of SL data, enabling us to measure the 
dark matter mass profile shapes and mass concentrations 
from hundreds of multiply-imaged sources, providing pre- 
cise (~ 10%) observational challenges to scenarios for the 
DM mass di stribution (for full det ails about the CLASH 
program see jPostman et aL]|2011| . 

The 16 HST bands chosen tor this project ranging from 
the UV through the optical and to the IR, and additional 
spectra available from large ground-based telescopes for 
some of the brighter arcs, enable us to obtain accurate 
redshifts for the multiply-lensed sources presented in this 
work. We use these remarkable imaging data alo ng with 



our well-tested approach to SL modeling (e.g ., Broad 



hurst et al.|[2005b] |Zitrin et al.||2009a|b[ |2010| |20TTajc| 



in order to find a significant number of multiple images 
across the central field of Abell 383 (A383 hereafter) so 
that its mass distribution and profile can be constrained 
with high precision. Various other ma ss models for this 
cluster were previously presented (e.g.,|Smith et al. 2001 

,^ s \ ,~> — I 1 1- 1 -i ! l Innn jl Urt^^l In y ' . I I J\ f\ A -I I , ' II 



2005||Sand eTal . 2004, 2008 ; [Newman et al.|201l[ ) usually 



based on WFPC2/HST single-band observations, uncov- 
ering 3-4 multiple image-systems and various candidates, 
as will be further discussed in §4.1| 

The approach to SL modeling implemented here in- 
volves only six free parameters so that in practice the 
number of multiple images uncovered readily exceeds the 
number of free parameters as minimally required in or- 
der to obtain a reliable fit, allowing for identification 
of other multiply-lensed systems across the cluster field. 
Our approach to lens-modeling is based on the reason- 
able assumption that mass approximately traces light. 
We h ave independently tested this assumption in Abell 
1703 (|Zitrin et al.|2010P, by applying th e non-param etric 
technique of |Liesenborgs et al.| f2006| |2007[ |2009[ ) for 
comparison, yielding similar results. Sucn parameter- 
free methods usually do not have the precision to actu- 
ally find new multiple- images, but the resulting ID radial 
profiles are sufficiently accurate for meaningful compar- 
isons. Independently, it has been found that SL methods 
based on parametric modeling are accurate at the level 



CLASH Strong-Lensing Analysis of A383 



3 



of a few percent in determ ining the projected inner mass 
( |Meneghetti et al.||2010b| ). 

The paper is organized as follows: In §2 we describe the 
observations, and in §3 we detail the SL analysis. In §4 
we report and discuss the results where in §5 we compare 
these to numerical simulations. The results are then sum- 
marized in §6. Throughout this paper we adopt a con- 
cordance ACDM cosmology with (O mu = 0.3, ^ao = 0.7, 
h = 0.7). With these parameters one arcsecond corre- 
sponds to a physical scale o f 3.17 kpc for this cluster (at 
z = 0.189; |Sand et a l. 2004). The reference center of our 
analysis is fixed on the brightest cluster galaxy (BCG): 
RA = 02:48:03.41 Dec = -03:31:44.91 (J2000.0). 

2. OBSERVATIONS AND REDSHIFTS 

As part of the CLASH program (see §1), Abell 383 
was observed with HST between 2010 November to 2011 
March. This is our first of 25 clusters to be observed 
to a depth of 20 HST orbits in 16 filters with the Wide 
Field Camera 3 (WFC3) UVIS and IR cameras, and the 
Advanced Camera for Surveys (ACS) WFC. Observation 
details and filters are provided in Table [l] 

The images are processed for debias, nats, superflats, 
and darks, using standard techniques. The ACS im- 



TABLE 1 

CLASH HST Observation Log for A383 



ages are further corrected for bias striping (Grog in et al. 
20101) and C TE/CTI degradation effects QAnderson & 
Bedin| 12010) . WFC3/IR pixels are flagged and down- 
weighted tor persistence effects. All images are then co- 
aligned and combined using drizzle algorithms to a scale 
of 0.065"/ pixel. An additional set of images with the 
original ACS 0.05"/ pixel scale is produced, onto which 
we apply our modeling initially to maintain the higher 
resolution, where the full UVIS/ ACS/ WFC3-IR data set 
is then importantly used for multiple-images verification 
and measurement of their photometric rcdshifts. Further 
details of our pipeline will be presented in an upcoming 
paper. 

Based on the 16-filter photo metry, we obtain pho- 
tometric redshifts using BP Z (Bemtez 2000 Bem'tez 
et al .||2004| |Coe et al.||2006D and LePhare (LPZ here- 
after 



Arnouts et al 
two methods yieldec 



1999 



some 



llbert et 
ot the 



al 
best 



2006). These 
results of all 



photo-z m ethods tested by the PHo to-z Accuracy Test- 
ing group ( |Hildebrandt et al.||2010) . BPZ and LPZ are 
similar in that spectral energy distribution (SED) tem- 
plates are redshifted and fit to observed photo metry. 
BPZ currently uses 6 tem plates from PEGASE ( |Fioc| 
& Rocca-Volmerange 19971, calibrated using the FIRE 
W(JRKS~p hotometry and spectroscopic redshi fts from 
Wuyts et al.| (|2 008). LPZ uses templates fr om Bruzual 
&; Cha riot ( 2003]rcahbrated using COSMOS ( |Koekemoer 



et al. |2UU7 ) photometry a n d spe ctroscopic redshifts as 
described 



llbert et al. 



(|2009j). The templates are 
empirically generated, describing well the full range of 
galaxy colors found in these multiband catalogs (less 
than ~ 1% outliers for high quality spectroscopic sam- 
ples), and therefore implicitly encompass all the range of 
metallicities, extinctions, and star formation histories of 
real galaxies. Further details on these methods can be 
found in the aforementioned references. While similar, 
the two methods serve as important cross-checks of one 
another. 

The distances to the galaxies are, of course, key in- 
gredients to the lens model. The photo-z analyses used 
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Note. — Observation were carried out between 2010, November 
18th, to 2011, March 3rd. The table summarizes the total exposure 
time in each filter. Note that these values are specific to A383. 
Observation times may vary for other CLASH clusters. We also 
note that the 5<r limiting magnitude is fainter than 26.8 AB mag 
for all 16 filters, as will be detailed in an upcoming paper ^Postman) 
|et al.|2011| . 



here also clearly aid us in assessing the robustness of 
the multiple-image identifications. The photometry of 
some lensed images may be significantly contaminated 
by brighter nearby cluster galaxies. SExtractor attempts 
to correct for this by measuring and subtracting the lo- 
cal background around each object. This works well in 
some cases but not all. To better reveal these lensed 
galaxies, we have carefully modeled and subtracted the 
light of several cluster galaxies including the BCG. While 
this improved the detection of some lensed galaxies, it 
did not consistently improve their photometry and thus 
photometric redshifts. The cluster galaxy wings must be 
modeled and subtracted very robustly and consistently 
to achieve quality photometry in all 16 bands for faint, 
nearby galaxy images. 

Explicitly, for its subtraction , the BCG has been mod- 
eled using the CHEF basis ( Jimenez- Teja & Bemtez 
This basis comprises both Chebyshev rational and 



2011) 



trigonometric functions, ensuring that the extended disk 
of this object is properly modeled. The flexibility of the 
CHEFs scale parameter allows us to accurately represent 
the BCG while keeping significant substructure and arcs 
unchanged. An example of the BCG subtraction is seen 
in Figure [3] 

Our 16 niters were selected based on tests with sim- 
ulated ph otometry to yield pre cise (Az ~ 0.02(1 + z)) 
photo-z 's ( |Postman et al.|[2011[ ). Previous work has also 
demonstrated how photo-z precision improves by increas- 
ing the number of (prefe rably overlapping) fi lters for a 



fixed total observing time pBemtez et al. 2009b ) . The em- 
pirical precision of CLASH photo-z's for arcs and other 
galaxies, including the relative contributions of various 
filters, will be detailed in future work. 

3. STRONG LENSING MODELING AND ANALYSIS 

We apply our well tested approach to lens model- 
ing, which has previously uncovered large numbers of 
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1000 1500 
ACS pixels 



Fig. 1. — The starting point of the mass model, where we define 
the surface mass distribution based on the cluster member galaxies 
(see 83). Axes are in ACS pixels (0.05" /pixel), and a 20" bar is 
overplotted. North is up, East is left. 



multiply-lensed galaxies in ACS images of Abell 1689, 
C10024, 12 high-z MACS clusters, MS 13 58, and the 
"Pandora cluster" Abell 2744 (respective l y, [B roadhurst 
et~aT1|2005bl |Zitrin et al.1|2009b[ |2011a|c[ fMerten et aT 
201ip . Briefly, the basic assumption adopted is that mass 
approximately traces light, so that the photometry of the 
red cluster member galaxies is used as the starting point 
for our model. Cluster member galaxies are identified 
as lying close to the cluster sequence by the photometry 
described in £[2j In addition, using our extensive 16-band 
imaging and corresponding photometric redshifts, these 
can be then verified as members lying at the cluster's 
redshift. 

We approximate the large scale distribution of clus- 
ter mass by assigning a power-law mass profile to each 
galaxy (see Figure[T]), the sum of which is then smoothed 
(see Figure [2]). The degree of smoothing (S) and the in- 
dex of the power-law (g) are the most important free 
parameters determining the mass profile. A worthwhile 
improvement in fitting the location of the lensed images 
is generally found by expanding to first order the gravi- 
tational potential of this smooth component, equivalent 
to a coherent shear describing the overall matter ellip- 
ticity. The direction of the shear (0 7 ) and its amplitude 
(I7I) are free parameters, allowing for some flexibility in 
the relation between the distribution of DM and the dis- 
tribution of galaxies, which cannot be expected to trace 
each other in detail. The total deflection field cEx(#), 
consists of the galaxy component, a ga i(6), scaled by a 
factor Kg a i, the cluster DM component aDM{9), scaled 

by (1-Kg a i), and the external shear component a ex {6), 
all scaled by the overall normalization factor K q : 



Fig. 2. — The resulting smooth mass component of the mass 
model (see §311. Axes are in ACS pixels (0.05" /pixel), and a 20" 
bar is overplotted. North is up, East is left. 

where the deflection field at position 6 m due to the ex- 
ternal shear, a ex {8 m ) = (a eXtX ,a ex>y ), is given by: 

) = |7| cos(20 7 )Aa; m + I7I sin(2</> 7 )Ay m , (2) 



a T (0) = K q {K gal a ga i{6) + (1 - K gal )a DM (9) + - 



) = I7I sin(20 7 )Aa; m - |-y | cos(2</> 7 )Ay m , (3) 

where (A%, Ay m ) is the displacement vector of the po- 
sition 9 m with respect to a fiducial reference position, 
which we take as the lower-left pixel position (1, 1), and 
7 is the position angle of the spin-2 external gravita- 
tional shear, measured counter-clockwise from the cc-axis. 
The normalization of the model (K q ) and the relative 
scaling of the smooth DM component versus the galaxy 
contribution (K ga i) bring the to tal number of free pa- 
rameters in the model to 6 (see Zitrin et al. 2009b for 
more details). This approach to SL is sufficient to ac- 
curately predict the locations and internal structure of 
multiple images, since in practice the number of multi- 
ple images uncovered readily exceeds the number of free 
parameters, so that the fit is fully constrained. 

In addition, two of the six free parameters, namely the 
galaxy power law index q, and the smoothing degree <!?, 
can be primarily set to reasonable values so that only 4 
of the free parameters have to be constrained initially, 
which sets a very reliable starting-point using obvious or 
known systems. This is because these two parameters 
control the mass slope, but the overall mass distribu- 
tion and corresponding critical curves do not strongly 
depend on them. The mass distribution is therefore 
primarily well constrained, uncovering many multiple- 
images which can then be iteratively incorporated into 
the model, by using their redshift estimation and loca- 
tion in the image-plane. 

We use this preliminary model to delens the more obvi- 
ous lensed galaxies back to the source plane by subtract- 
ing the derived deflection field. We then relens the source 
o(0)), plane in order to predict the detailed appearance and lo- 
(1) cation of additional counter images, which may then be 
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Fig. 3.— Galaxy cluster A383 (z = 0.189) imaged with HST/ACS/WFC3. North is up, East is left. We number the multiply-lensed 
images used and uncovered in this work. The numbers indicate the 27 lensed images, 13 of which correspond to 4 newly identified sources, 
and the different c olors are used to distinguish the 9 different sources. For more details on the each system and the robustness of the new 
identifications see §4.1| The overlaid white critical curve corresponds to systems 3 and 4, at z s = 2.55, enclosing a critical area of an effective 
Einstein radius of ~ 52 kpc at the rcdshift of this cluster (16.3"). Also plotted is a red critical curve, which corresponds to system 5, 
the drop-out high redshift galaxy at z s = 6.027. The composition of this color image is Red=F105W+F110W+F125W+F140W+F160W, 
Green=F606W+F625W+F775W+F814W+F850LP, and Blue=F435W+F475W. This image was generated automatically by using the 
freely available Trilog y softwareF] The upper-right inset shows the central core with the BCG subtracted, using the method of |Jimenez-| 
| Teja fc Bem'tez| l |201l| see also q2T. 
3 http: //www.stsci.edu/~dcoe/trilogy/ 



identified in the data by morphology, internal structure 
and color. The best fit is assessed by the minimum x 2 
uncertainty in the image plane: 

X 2 ^Y,^ x -x t f + {y t - yi f)/a\ (4) 

i 

where x i and y i are the locations given by the model, 
Xi and yi are the real image locations, a is the error in 
the location measurement (taken as 0.5"), and the sum 



is over all N images. The model location of each image 
is the averaged location given by relensing all other im- 
ages of the same system. The best-fit solution is unique 
in this context, and the model uncertainty is determined 
by the location (of predicted images) in the image-plane 
itself. Importantly, this image-plane minimization does 
not suffer from the bias involved with source-plane min- 
imization, where solutions are biased by minimal scat- 
ter towards shallow mass profiles with correspondingly 
higher magnification. 
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The model is successively refined as additional sets of 
multiple images are incorporated to improve the fit, im- 
portantly using also their redshift information for bet- 
ter constraining the mass slope. The mass profile is 
coupled to the redshift distribution of the different sys- 
tems, since for each redshift the enclosed mass and cor- 
respondingly the deflection angle, depend on the lens 
and source angular-diameter distances (£>;, D s , respec- 
tively). Explicitly, the deflection angle is defined as 

ot(0) = 4G ^,2g <e ' > j 1 ^, and since the lens distance is con- 
stant, the mass slope is constrained through the cosmo- 
logical relation of the Di s /D s growth with source red- 
shift, where Di s is the distance between the lens and the 
source. This is seen more clearly in Figure [T2} 

4. RESULTS AND DISCUSSION 
4.1. Multiple- Images, Mass Model and Critical Curves 



In addition to t he previously-known sys tems (see New- 
Richard et al. 



man et aL] 2011 



2011) 



and references therein 
our modeling technique has uncovered 13 new 
multiply-lensed images and candidates in the central field 
of A383, belonging to 4 new systems. We thus substan- 
tially increase in this work the number of available con- 
straints on the mass profile of this cluster. 

We have made use of the location and redshift infor- 
mation of the multiple-images to fully constrain the mass 
model. In our minimization procedure, we obtain for 
most important parameters controlling the mass distri- 
bution, values of q = 1.08 ± 0.08 and S = 12 ± 2, but 
note these are highly coupled to the photometry used to 
construct the mass model, and to our procedure detailed 
in §3. 

We find that the critical curve for a source at z s = 2.55 
(systems 3-4) encloses an area with an effective Einstein 
radius of r^ = 16.3 ± 2", or ~52 kpc at the redshift 
of the cluster. A projected mass of M — 2.4 ± 0.2 x 
10 13 M Q is enclosed by this critical curve (see Figure 
For general comparison, this is in good agreement wit 
the Einstein rad ius-mass relation f or a source at z s ~ 
2 — 2.5, found in Zitrin et al.||2011a| (taking into account 



a- 

ith 



also the different lens distances; see Figure 27 therein). 
This is naturally expected from the lensing equations, 
though constitutes an important consistency check. The 
corresponding critical curves are plotted on the cluster 
image in Figure|3]along with the multiply-lensed systems. 
The resulting mass distribution and its profile are shown 
in Figures H] and [5] 

It shoula be stressed that the multiple-images found 
here are accurately reproduced by our model and are not 
simple identification s by eye. The parametric method of 



Zitrin et al. 



([2009b) has been shown in many cases to 
have the predictive power to find multiple images in clus- 
ters. Due to the small number of parameters this model 
is initially well-constrained enabling a reliable identifi- 
cation of other multiple-images in the field, which can 
be then used to fine-tune the mass model. Naturally, 
the mass model predictions have to be identified in the 
data and verified further by comparing the SEDs and 
photometric redshifts of the candidate multiple-images, 
especially in cases where the images are not prominently 
bright and big, so that internal details cannot be reliably 
distinguished. As some of the objects identified here are 
faint and some may be contaminated by nearby cluster 



members even after their subtraction, for the less secure 
cases we supply also the photo- z distributions and spec- 
tral energy distributions (SEDs) from our 16 HST-band 
imaging, so that the reader could assess the plausibility of 
these identifications. We now detail each multiply-lensed 
system, as listed in Table [2] 

Systems 1 — 2 : The prominent giant arc, nearly 20" 
long, most likely consists of two sources (systems 1 and 2 
here) at the same r edshif t of z s = 1.01 QSand et al.|2004 



20081 Smith et al. 2005). An additional radial counter 



image is seen m the UCG halo (see also [ Newman et al 
2011 1. This system was identified by |Smffn~eTal.| ( |200l| ) 



in WFPC2 1-band imaging, who also spectroscopically 
measured the west side of the arc to be at z s — 1.01. 



Following measurement of Sand et al. (2004) with a slit 
passing through the BCG, the radial arc, and the eastern 
part of the main arc, yielded an identical redshift of z s = 
1.01 for both as well. 

Following examples from other well known clusters, it 
is not common for a giant arc to consist of two different 
sources. We therefore primarily do not use the location of 
the multiple- images of these systems in our minimization 
(only their redshift), though our model agrees with this 
previous interpretation and accurately produces these 
multiple-images at this redshift. In addition, our model 
suggests that part of the radial arc is also contributed 
by the left side of the giant arc (system 2). However, it 
is still plausible that the giant arc consists of only one 
elongated source. We find that the full giant arc, when 
projected back to the source plane, corresponds to ~11 
kpc in length, which may indeed be accounted for by a 
single source. The reproduction of the source is seen in 
Figure [6[ 

Systems 3 — 4 : Corresponding to a pair of sources, 
the images of which are lensed to appear next to each 
other on the two sides of a prominent cluster galaxy, as 
can be seen in Fig ure [3[ These images were identified by 



Smith et al. ( 2001 ) who mapped the internal str uctures in 
detail , and w ere spectroscopically measured by [Newman] 



et al. ([20111) to be at a redshift of z s = 2.55. In addition, 
our IK/ WFC3 images show clearly, for the first time, that 
these two systems indeed have two different colors and 
SEDs, and our mass model accurately reproduces each as 
shown in Figure [7] We note that our model suggests that 
a small part of the radial arc may consist of a counter 
image of these systems, in addition to the radial images 
of systems 1 and 2. 

Eastwards to image 3.3 there is a faint arc which might 
be related either to this system or to system 4. This faint 



extending ar c was marked as part of this system by Smith 



et al. 



2001) but omitted in recent analysis (Newman 



et al.|201l| ). We find that this faint arc, marked as 3.5/4.4 



here, may be related to this system (also yielding a sim- 
ilar photometric redshift of ~ 2.7), and is reproduced as 
part of this system by our model if we slightly increase 
the weight of its neighboring galaxy (RA=02:48:03.42, 
DEC=-03:32:09.02), see Figure^} On the other hand, the 
IR colors do not strongly support connection to systems 
3 and 4, and this faint arc might be a locally multiply- 
lensed separate system. In any case its inclusion has only 
a negligible and local effect on the mass model. 

System 5 : Two images of a multiply-lensed Lyman- 
break, hi gh redshift galaxy at z s — 6.027, reported re- 
cently by [Richard et al. (2011) based on CLASH imag- 



CLASH Strong-Lensing Analysis of A383 



7 



TABLE 2 
Multiple-image systems 



ARC 




RA 


DEC 


RP7 7 

^phot 


T P7 7 


spec-z 


■^model 


Comment 


ID 


(J2000.0) 


(J2000.0) 


(best) [95% C.L.] 


(best) [99% C.L.] 








1.1 


no 




no qqi 


-03:31:49.72 


0.97 [0.83-1.05] 


0.93 [0.92-0.93] 


1.01 


(I oil 




1.2 


09 

UZ 


TO 


uo.ozo 


-03:31:41.85 


0.53 [0.34-0.59] 


0.47 [0.47-0.48] 


1.01 




radial image in BCG halo 


2.1 


no 


-lo 


no nyi7 


-03:31:58.95 


0.95 [0.87-1.03] 


0.90 [0.90-0.92] 


1.01 


d oil 




2.2 


02 


48 


02.852 


-03:31:58.04 


0.96 [0.82-1.04] 


0.85 [0.78-0.93] 


1.01 






2.3 


no 


-1 Q 

■lo 


no a c;o 


-03:31:52.84 


0.84 [0.77-0.91] 


0.76 [0.67-0.84] 


(1.01) 


„ 




3.1 


02 


■48 


02.426 


-03:31:59.40 


2.79 [2.64-2.94] 


2.90 [2.54-3.15] 


2.55 


(2.55) 




3.2 


02 


■48 


no onn 
02.309 


-03:31:59.21 


2.90 [2.75-3.05] 


3.01 [2.92-3.08] 


2.55 






3.3 


no 


lo 


nQ noft 


-03:32:06.75 


2.56 [2.42-2.70] 


3.03 [2.86-3.09] 


2.55 


■■ 




3.4 


no 


lo 


no Qnn 


-03:32:01.74 


2.88 [2.73-3.05] 


3.01 [2.87-3.16] 


(2.55) 


„ 




4.1 


02 


18 


no O A A 
02.244 


-03:32:02.07 


0.20 [0.15-0.25] 


0.20 [0.20-0.26] 


2.55 


(2.55) 




4.2 


no 


A S2 
10 


no 01 a 

UZ.Z14 


-03:32:00.25 


2.85 [2.70-3.00] 


2.91 [2.82-3.01] 


2.55 






4.3 


no 


-1 Q 
10 


no Q/1 7 


-03:32:06.68 


3.09 [2.93-3.25] 


3.05 [2.91-3.20] 


2.55 


„ 




5.1 


no 


10 


ni ok/1 


-03:31:34.77 


5.95 [5.68-6.22] 


5.87 [5.64-5.99] 


6.027 


(6.027) 




5.2 


no 


a a 
lo 


ha find 
U4.0UU 


-03:31:58.47 


6.01 [5.74-6.29] 


5.96 [5.72-6.12] 


6.027 






6.1 


02 


48 


04.272 


-03:31:52.77 


2.67 [2.53-2.81] 


2.13 [2.05-2.18] 




~ 2.0 




6.2 


02 


48 


03.377 


-03:31:59.27 


2.38 [2.25-2.51] 


1.93 [1.90-2.09] 








6.3 


02 


48 


02.153 


-03:31:40.88 


1.89 [1.78-2.04] 


2.10 [1.90-2.20] 


_ 






6.4 


02 


48 


03.720 


-03:31:35.87 


1.80 [1.69-1.91] 


1.54 [1.46-1.57] 






bright galaxy nearby 


7.1 


02 


48 


04.089 


-03:31:25.54 


4.60 [0.64-4.82] 


4.50 [4.24-4.76] 




~ 4.6 


bimodal 


7.2 


02 


48 


03.568 


-03:31:22.55 


4.65 [0.38-5.20] 


4.77 [0.52-5.58] 








7.3 


02 


48 


03.130 


-03:31:22.16 


4.70 [4.35-5.07] 


4.56 [0.20-5.08] 








8.1 


02 


48 


03.681 


-03:31:24.43 


0.34 [0.24-2.43] 


0.33 [0.20-3.19] 




~ 3.1 


bimodal 


8.2 


02 


48 


03.386 


-03:31:23.46 


2.94 [2.41-3.25] 


2.93 [0.20-3.46] 








9.1 


02 


48 


03.920 


-03:32:00.83 


3.91 [3.63-4.10] 


3.83 [3.51-4.09] 




~ 4.0 




9.2 


02 


48 


04.046 


-03:31:59.21 


0.48 [0.26-0.54] 


0.47 [0.39-0.53] 






segment yields z„h ot ~ 3.9, see f 


1.1 




9.3 


02 


48 


03.872 


-03:31:35.03 


3.96 [3.77-4.15] 


3.57 [3.56-3.59] 








9.4 


02 


48 


01.918 


-03:31:40.23 


3.80 [3.61-3.99] 


3.75 [3.57-3.82] 









Note. — Multiple-image systems and candidates used and uncovered by our model. For more detailed information on each system see 
the corresponding subsection. The columns are: arc ID; RA and DEC in J2000.0; best photo-z using BPZ, along with 95% confidence level, 
minimal and maximal photo-z; best photo- ,2 using LPZ, along with 99% confidence level, minimal and maximal photo-z; spectroscopic 
redshift, spec -z ; z^nj^i, e s timated red shift for the arcs which lack spectroscopy as predicted by the mass model; comments. System 1 was 
uncovered by |Smith et al.| |2001 2005) who measured its redshift spectroscopicallv , which is the value given below. Systems 2-5 were also 
found and spectroscopically measured in previous works (Smith et aL||200l| |2005| [Sand et al^|2004| |2008| |Newman et al ||2011| |Richard| 
|et al.|2011| ). Note also that unusually large errors in the photo-z imply a bimodal distribution. In such cases the values which agree with 
the bL model can be different than specified in the best photo-z column, as they arise from another peak in the distribution. Such cases 
are specified in the comments. 



ing and Keck spectra. We also identify these two images 
and measure photometric redshifts of z ss 6.01 and 5.95. 
The high redshift of this system expands substantially 
the lensing-distance range thus enabling us to c onstr ain 
the profile with better accuracy, as discussed in §4.2| 

System 6 : This system consists of 4 blue images with 
similar internal details including a brighter white blob, at 
a typical photometric redshift of z s ~ 2.4 for this system 
(see Table M Fi gure [9L Our model reproduces these im- 
ages very well (Figure 8 ) , though it slightly favors a lower 
redshift of z s ~ 2 but aue to the distances involved this 
is in practice only a ~ 1% difference in the redshift dis- 
tance ratio. These images were matched up for the first 
time in this work enabled by the deep, high- resolution 
HST data. Due to the variance in the SEDs and there- 
fore photometric redshifts of the images of this system, 
we supply also the photo-z distributions and SEDs in 
Figure [9j so that the reader could more easily assess the 
plausibility of this system. 

Systems 7 — 8 : Two thin and long arcs following simi- 
lar symmetry, at a relatively high redshift of z ~ 4.5 and 
2^3, respectively. Their symmetry especially with re- 
gards to the critical curves, shows beyond a doubt that 
these are multiply-lensed systems (see also Figure 10), 



despite being too faint to measure their photometric red- 
shift unambiguously. These images as well were matched 
up for the first time in this work. 

System 9 : A faint, wide greenish- looking arc 17" south 
east of the BCG (see Figure [3]). Our model accurately 
reproduces this arc as a double image. In addition, two 
other small counter-images are predicted, for which we 
identify the best-matching candidates in the data. These 
images were matched up for the first time in this work, 
and except for image 9.2, sho w s imilar photometric red- 
shifts of ~ 3.8 (see also Figure 11 ), in agreement with our 
model prediction. In addition, it should be noted that 
photo-z analyses of some segments of the arc designated 
as 9.2 imply indeed a redshift of ~ 3.8, similar to the 
other three images of this system. We also acknowledge 
the possibility that other similar looking objects near-by 
images 9.3 and 9.4 may be the actual counter images - 
especially since 9.3 and 9.4 seem slightly brighter than 
9.1 and 9.2. Such a degeneracy however does not affect 
the mass model in a noticeable way. Due to the variance 
in the SEDs and therefore photometric redshifts of these 
images, we sup ply also their photo-z distributions and 
SEDs in Figure [Tl] so that the reader could more easily 
assess the plausibility of this system. 
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Fig. 4. — 2D surface mass distribution (re), in units of the critical 
density (for z a = 2.55), of A383. Contours are shown in linear 
units and in spaces of Are = 0.1, derived from our mass model 
constrained using the many sets of multiply-lensed images seen 
in Figurej3] As can be seen, the mass distribution is fairly round. 
Axes are in ACS pixels (0.05" /pixel), and a 20" bar is overplotted. 
North is up, East is left. 
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Fig. 5. — Radial surface mass density (re) profile in units of 
the critical surface density (for a source redshift of z s = 2.55). 
The inner SL data were derived using the sets of multiple im- 
ages shown in Figure [3] We overplot our preliminary WL data 
analysis. As can be seen, these are in very good agreement in 
the region of overlap. A joint SL+WL gNFW fit yields M vir = 
(5.94to;g7 ± 0.71) x W 14 M Q /h (or M vir ~ 8.49 x 1O 14 M ) and a 
concentration parameter of c_2 = 7.95^g'gp ± 0.55. The parame- 
ter a = 1.078^q'qy3 0-059, so that the overall fit is similar to a 



simple NFW (see j ]4.2| for explicit comparison) . These values are 
in common with more massive well studied cl uste rs, and lie above 
the standard c-M rel ation, as seen in Fi g. |13| Also plotted is 
the ID WL analysis of |Huang et al.|(j2011^ . A clear consistency is 
seen through the extensive W-L range, tnough our profile is more 
consistent with the SL data and is not underestimated in the inner 
region. A more thorough, 2D WL analysis will be published soon 
(Umetsu et al., in preparation). 

4.2. Mass Profile 

The inner mass profile is accurately constrained by in- 
corporating the cosmological redshift-distance relation, 
i.e., the lensing distance of each system based on the 
measured spectroscopic or photometric redshifts. In so 
doing we normalize our mass model to systems 3 and 4, 




Fig. 6. — Reproduction of systems 1 and 2 (consisting the giant 
arc) by our model. We lens the full giant arc to the source plane and 
back, with a lensing distance corresponding to z s = 1.01. Over- 
plotted are the reproduced radial arc, 1.2, and an enlarged image of 
the reproduced source along with its physical scale. The prominent 
bright blob lensed in the procedure may be unrelated. The source 
image had color manipulation and noise cleaning procedures acted 
on, to better show the internal details. 





Fig. 7. — Reproduction of systems 3 and 4 by our model, by 
delensing jointly images 3.1 and 4.1 into the source plane with a 
lensing distance corresponding to z s = 2.55, and then relensing 
the resulting source plane pixels onto the image plane. Note the 
different colors of these systems, seen in an WFC3/IR color image. 

so that the normalized scaling factor, f(di s /d s ), is equal 
to 1 for z s = 2.55. We then make use of the z = 1.01 
system, and the highest-z system at z s = 6.027, in or- 
der to expand the f{di s /d s ) range, along with the other 
systems whose photometric or spectroscopic redshifts are 
incorporated to constrain the profile. The resulting mass 
profile is seen in Figure 15] 

We examine how wellthe cosmological relation is re- 
produced by our model, accounting for all systems with 
spectroscopic or photometric redshifts, as shown in Fig- 
The predicted deflection of the best fitting model 



12 



at the redshift of each of these systems clearly lies along 
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Fig. 8. — Reproduction of system 6 by our model, by delensing 
image 6.2 into the source plane with a lensing distance correspond- 
ing to z a ~ 2, and then relensing the resulting source plane pixels 
onto the image plane. Our model clearly reproduces accurately the 
other images in this system. 
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Fig. 9.— Photo-z distribution (blue) and 16 HST-band SEDs 
of the four images of system 6 (see Figure |8]l , generated using 
BPZ (Q. As can be seen, the photo-z's support in general our 
identification of this system, though with some uncertainty and 
corresponding variation in the SEDs. The yellow stripe in the 
photo-z distributions corresponds to the cluster redshift at z = 
0.189, and the red crosses in the SED plots mark the lcr magnitude 
detection limit in each filter. 



Fig. 10. — Stamp images of some of our newly identified multiple- 
images, seen more clearly in this color-composite image constructed 
from 16 HST bands ranging from the UV to the IR. As can be 
seen, systems 7 and 8 follow the same symmetry and are clearly 
multiple images as supported also by our mass model and photo-z's 
(see Table|2l. In system 9, images 9.1 and 9.2 are clearly multiply- 
lensed to form the greenish-looking arc, while 9.3 and 9.4 are the 
most likely counter ima ges a s predicted by our model and coherent 
photo-z's seen in Figurc flT] Other similar looking objects are seen 
close to these images which might be the actual counter images, but 
such a degeneracy does not affect the mass model in a noticeable 
way. 

the expected cosmological relation, wi th a small mean 
deviation of only A/ < 0.01 (see Figure 12), strengthen- 
ing the determination of the mass profile slope. 

In addition, we note that our mass profile shows con- 
sistency with a r ecent joint lensing, X- ray ,and kine- 



matic analysis by Newman et al. ( 201lj as read from 
Fig. 2 therein), out to at least twice the Einstein ra- 
dius where our SL data apply For example, for the ra- 
dius of the giant tangential arc (systems 1 and 2), the 
model of Newman et al. encloses a projected mass of 
~ 2 x 1O 13 M0, while our model yields for that radius 
~ 2.2 x 1O 13 M . At higher radii, say a 100 kpc (which is 
about twice the Einstein radius), both models yield simi- 
larly ~ 6x 1O 13 M . Due to the different interpretation of 
the radial arc, some differences are seen in the very inner 
region, so that for radii of 5-10 kpc (1.5 — 3") our model 
yields - 0.09-0.2 xl0 13 Af Q , versus - 0.05-0. 7x 10 13 M Q 
for the model by Newman et al. 

We combine our SL-based profile with ID WL distor- 
tion and magnification measurements out to and beyond 
the virial radius (R V i r — 11.3 arcmin; or ~2.1 Mpc; cor- 
responding to an overdensity of ~115 with respect to the 
critical density of the universe at the cluster redshift), 
obtained from deep multicolor Subaru imaging (see Fig- 
ure [5]). Here we have chosen the BCG position as the 
center of mass for our mass profile analysis, where our 
strong-lens modeling shows that the dark-matter center 
of mass is consistent with the location of the BCG, with- 
out any noticeable offset within errors. The SL profile 
is obtained in 81 linearly-spaced radial bins from 9 = 2" 
(excluding the BCG) to 42", including cosmic covariance 
between radial bins due to the uncorrelated large scale 
structure, estimated by projecting the nonlinear matter 
power spectrum out to the median depth of z s = 2.55 (see 
Table 2), following the prescription detailed in Umetsu 



et al.| ( |2011a l 

The WL mass profile, given in logarithmically-spaced 
radial bins, was derived using the Bayesian method of 
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Fig. 11. — Another example of the photo-z distribution (blue) 
and 16 HST-band SEDs, generated using BPZ and used to 

match the four images of system 9 (see Figure [3^ The photo-z 
distributions and the SEDs support the identification of this system 
at z s ~ 3.8, apart from image 9.2, although note that photo-z 
analyses of some of its segments do favor a higher-z of z s ~ 3.8, 
similar to the other images in this system. The yellow stripe in 
the photo-z distributions corresponds to the cluster redshift at z = 
0.189, and the red crosses in the SED plots mark the lc magnitude 
detection limit in each filter. 



Umetsu et al. ( 2011a|b l that combines WL tangential- 
distortion and magnification-bias measurements in a 
model-independent manner, with the assumption of 
quasi-circular symmetry in the projected mass distri- 
bution]^ The method applies to the full radius range 
outside the Einstein radius, and is free from the mass- 
sheet degeneracy, recovering the absolute mass normal- 
ization or equivalently the projected mass M 2 d{< #min) 
(corresponding to the first WL bin of Figure 5) inte- 
rior to the inner radial boundary of WL measurements, 
(?min = 42" r^). The strong and weak lensing are in 
excellent agreement where the data overlap, 8 ~ 0.7— 1' 
85 - 190 kpc). 
For co mparison we oyerplot in Figure [5] also the recent 



profile of Huang et al 
WL distortion data. T 



(2011) derived from the Subaru 



re two profiles are in good agree- 



24 This method applies without the axial symmetry approxima- 
tion in the WL regime where nonlinearity between the surface mass 
density and observables is negligible. 



ment and very similar in the WL regime, but the Huang 
et al. profile is slightly underestimated in the inner region 
relative to our SL data. Our secure ba ckground selection 



method ( Medezinski et aLpOlO 2011 ) carefully combines 
all color and clustering information to identify blue and 
red background galaxies in color-color space (B — R c vs. 
Rc ~ z'), minimizing contamination by unlensed cluster 
and foreground galaxies. It is important to stress that 
combining independent weak and strong lensing allows 
us to recover the full radial profile and ensure internal 
consistency in the region of overlap. 
We consider a gene ralized parametrization of the NFW 



( Navarro et al.||1996 ) m odel of the following form ( Zhao 



996| |Jing fc Suto 2000] ) 



PnfwM 



(r/r s ) a (l + r/r s 



(5) 



where p s is the characteristic density, r s is the character- 
istic scale radius, and a is the inner slope of the density 
profile. This model has an asymptotic outer slope of 
73£>( r ) = dhi p/dlnr = —3 (r — > oo), and reduces to the 
NFW model for a = 1. 

We refer to the profile given by equation [5] as the gen- 
eralized NFW (gNFW, hereafter) profile. It is useful to 
introduce the radius r_ 2 at which the logarithmic slope 
of the density is isothermal, i.e., 73D = —2. For the 
gNFW profile, r_ 2 = (2 — a)r s , and thus the correspond- 
ing concentration parameter reduces to c_2 = J"w/f*-2 = 
c V i r / (2 — a). We specify the gNFW model with the cen- 
tral cusp slope, a, the halo virial mass, M V i r , and the 
concentration, c_2 — c V i r /{2 — a). 

The joint SL+WL NFW fit yields M v „ = (5.37±g;^ ± 
0.26) x 10 14 M Q //i (or M v 
centration parameter of c V i r = 8.77 



7.67xl0 14 M m ) andacon- 



+0.44 

-0.42 



± 0.23, with a 



minimized \ 2 (Xmin) vaiue of 78.7/90 with respect to the 
degrees of freedom (dof), corresponding to a goodness- 
of-fit of Q — 0.798. Note that the values quoted include 
the statistical followed by the systematic uncertainty at 
a 68% confidence level. The systematic errors were esti- 
mated by changing the outer radial boundary of SL bins 
from 9 = 42"(= 6^) to 2r B (~ 33"). 

Only a very slight improvement in the fit is obtained 
by implementing the gNFW form described in eq. [5j 
A joint SL+WL gNFW fit yields M mr = (5.941^ ± 
0.71) x 10 14 M©//i (or M mr ~ 8.49 x 1O 14 M ), a con- 
centration parameter of c_2 = 7.951q'9q ± 0.86, and 

a = 1.078±g$| ± 0.059, with xLn/ dof = 77.5/89 and 
Q = 0.802. The central cusp slope a is consistent with 
unity, so that the overall fit is similar to a simple NFW, as 
also shown by the quoted x 2 an d Q values. These results 
are consistent with the values quoted by |Huang et al.| 
(M mr = 5.28t 4 ;^ x 10 14 M Q /h,c mr =5.68^), 



but the concentration is higher for example than the 
Chand ra X-ray based gNFW fit by |Schmidt fc Allen| 
( |2007| . 

We find that A383 lies above the standard c-M rela- 



tion (Figure 13 ), similar to several other well-known clus 
ters for which detailed lensing-based mass profiles have 
been constructed, addin g to the claimed tension with the 



standard ACDM model (Broadhurst et al.|2008|rUmetsu 
let al. 1120101 |2011b| |Zitrm et al.pOiOj see also ISadeh k 

Illivel 



Rephaeli 2008J). Still, the overall level of systematic un- 
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Fig. 12. — Growth of the scaling factor f(di a /d s ) as a function 
of redshift, normalized so f=l at z = 2.55. Plotted lines are the 
expected ratio from the chosen specified cosmological model. The 
circles correspond to the multiple-image systems reproduced by our 
mass model, versus their real spectroscopic or photometric redshift. 
The data clearly follow well the relation predicted by the standard 
cosmological model. As can be seen, based on only one cluster and 
due to the low number of multiple images, it may be still hard to 
discriminate between different cosmologies solely with these data. 
With a sample of 25 CLASH clusters however, where most clusters 
are expected to have larger critical area and thus larger numbers 
of multiple-images spread over a large redshift range, we will be 
able to put statistically-significant constraints on the cosmological 
parameters, combining the information from all clusters together. 

certainties may be too large to allow a definite conclusion 
regarding a clear inconsistency with ACDM predictions 
based on only a handful of clusters. This, in fact, is one 
of the primary goals of our CLASH program. Moreover, 
the WL data used here are based on ID analysis, while 
the concentration should be influenced by the triaxial- 
ity or other line-of-sight background structures, and this 
result will be revised in our following WL papers, using 
2D analysis (Umetsu et al., in preparation) and a joint 
SL+WL non-parametric reconstruction method (Merten 
et al., in preparation). 

In addition, we note that several dips are seen in our 
WL tangential distortion data in outer radii, resulting 
in positive perturbations in the k profile. We confirmed, 
by visual inspection in deep color Subaru imaging, that 
these correspond to several higher-redshif t background 
struc tures near the field of A383 (see also |Okabe et al. 
20101. In fact, the field is quite rich in such background 



structures, some of which are slightly magnified by the 
A383 foreground lens, as we will elaborate in our upcom- 
ing papers devoted to this configuration (Umetsu et al., 
Zitrin et al., in preparation). 

4.3. Brightest Cluster Galaxy 

Due to the presence of the radial arc in its halo, we 
may also constrain the mass enclosed within the BCG, 
and the corresponding M/L ratio. We find that the 
BCG encloses a projected mass of 1.14 ± 0.3 x 10 12 M Q 
within a radius of ~ 6" (~19 kpc) after subtracting the 
interpolated smooth DM component (~ 4.1 x 10 12 M Q 
inside this aperture). For comparison, recent stellar 



Fig. 13.— The joint SL+WL NFW fit of A383 (data point with 
lcr and 2cr confidence level contours) presented on the c—M plane, 
compared to expectations from simulations. Overplotted are the 
expected c-M rela tions and their lcr uncertainties, presented in 
|Duffy et al.| ( |2008| ) for their relaxed sample, scaled to z c = 0.189 
\blue band/. Also plotted are c —M relations for the full-sample 
clusters from|Duffy et al.| pOOS) ; |Zhao et aT] ( |2009^ |Klypin et 
(2010); Prada et al. (2011). As has been found tor other obscrvci 
clusters (e.g., Broadhurst et al. 2008 Oguri ct al. 2009 Zitrin et al. 
|2010f A383 appears to have a higher concentration than simulated 
clusters with similar masses and redshifts. We note that the WL 
data used here are based on ID analysis, while the concentration 
should be influenced by the triaxiality, and this result will be re- 
vised in our following WL paper using 2D analysis (Umetsu et al., 
in preparation). 

velocity-dispersion measurements of the BCG in A383 
yi eld a c± 450 ± 40 at t his radius (as read from Figure 2 
Newman et al.|[20TT |, which translates to a projected 



mass of 0. 891^15 x ^ va Mq , in agreement with our result 
We measure the BCG flux in several optical ACS 
bands, to obtain an average B-band luminosity of ~ 
8 ± 0.3 x 10 10 L o , within the aperture of ~ 6" (~19 kpc; 
fluxes were convert ed to luminosities usin g the LRG tem- 
plate described in Bemtez et al. 2009a). This yields a 
typical M/Lb of ~ 14 [M/L)q in this region, similar 
to other l ensing based BCG m asses in well-studie d clus 



ters (e.g., Gavazzi et aL]|2003| for MS 2137-2353, Zitrin 
fc Broadhurst||2009l for MACS J1149.5+2223, |Limousm 
et al.||2008[ |Zitrin et al.||2010| for Abell 1703), though of 



course the degeneracy between the BCG DM halo and 
the overall subtracted smooth cluster halo is still un- 
known. 

4.4. Modeling Accuracy and Uncertainty 

In general, since the deflection angle depends on the 
distance- redshift ratio (Di s /D s ), the SL modeling un- 
certainty, particularly with regards to the mass profile, 
is primarily coupled to the redshift measurement ac- 
curacy of the multiple-systems. In A383, five systems 
at three different redshifts have spectroscopic measure- 
ments, while the four other systems found in this work 
importantly supply four more constraints on the mass 
profile. Our 16-band ACS/WFC3 imaging, allows us to 
derive robust photo- z 's for all multiply- lensed systems 
discussed in this work, which we verify by using both the 
BPZ method, and the LPZ method (fjl. 

Still, due to the low number of parameters in 
our modeling, which constitutes a huge advantage for 
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finding multiple-images and producing efficiently well- 
constrained mass distributions, some local inaccuracy 
can be expected. Our best fit model for A383 reproduces 
all multiple-images described in this work within ~ 2.5" 
from their real location given their measured redshifts, 
aside from one image candidate belonging to system 9 
(see { 4.1 1 which is reproduced ~ 4" from its real loca- 
tion. In addition, we note that in our best-fit model, con- 
strained by all systems together, there is a slight offset 
of ~ 1" in the reproduced location of the radial arc, im- 
plying a slight inaccuracy in the BCG's very inner mass 
profile in that mode l. N ote however that this does not 
affect the result in |4.3| which was verified by comple- 
mentary models in which the radial arc is very accurately 
reproduced (but the fit is overall somewhat poorer taking 
into account all other systems), thus reliably constrain- 
ing the mass enclosed within the corresponding radius. 

The average image-plane reproduction uncertainty of 
our best-fit model is 1.68" per image in total, with an 
image-plane rms of 1.95" including all 27 multiply- lensed 
images. This image-plane rms is, for example, higher 
than that reported recently by Newman et al. (20111 
for A383 (rms = 0.27") based on only tour systems 
in two different spectroscopic redshifts, but is typical 
to most parametric-method reconstructions, when many 
mult iple-sy st ems are present. For example, Broadhurst] 
et al. ( 2005b ) a chieved an rms of 3.2" per image for A bell 
1689, and later Halkola et al. ( |2006 1 reported an rms of 
2.7" per image tor that cluster, while | Zitrin et al. ( 2009b I 



ant arcs (defined as having length-to- width ratios larger 



than 7 .5), were readily computed as in Meneghetti et al 
( 201 1 1 . To account for the lack of star formation in these 



produced an rms of ~ 2.5" for C10024. These values are 
comparable with our current model rms, taking into ac- 
count the difference in the critical area. 

In general, a higher number of parameters would sup- 
ply a more accurate solution, however the efficiency of 
a model and the confidence in it decrease substantially 
as more parameters are added to the minimization pro- 
cedure, especially if these are arbitrary non-physical pa- 
rameters as may be the case in other (non-parametric) 
methods. We have shown here as well as in many previ- 
ous examples (see also f|3| that our method, with a mini- 
mum number of fre e parameters, built on simple physical 
considerations (see Zitrin et al.||2009b] for full details), 
does a very good job in finding new multiply-lensed sys- 
tems, and thus in constraining the deflection field, and 
accordingly, the mass distribution and profile. 

5. COMPARISON TO NUMERICAL SIMULATIONS 

We now compare our derived A383 mass distribu- 
tion with cluster halos obtained from hydrodynamical 
simulations in the framework of the ACDM cosmology. 
The analysis we make h ere is inspired by the work of 
Meneghetti et al. (2011), where the Einstein ring sizes 
and the tensing cross sections of 12 massive MACS clus- 
ters modeled by Zitrin et al. (2011a) were compared with 
those expected from similar halos in the MareNostrum 



Universe cosmological simu l ation (Gottlober fc Yepes 
2007[ [Meneghetti et al.|2010a||Fedeh et al.|2010[ ). This is 
a 500 3 /i- J Mpc J volume filled with 1024 3 DM and 1024 3 
gas particles, evolved in the framework of a cosmological 
model with Q m = 0.3, — 0.7, and erg = 0.9. More 
details about th is simulation can be found in IGottloberl 
& Yepes (2007). For our comparison, we use halos ex- 
tracted from the same cosmological box, for which the 
median Einstein ring sizes and the cross sections for gi- 



simulations, we added to each halo a component mimick- 
ing the presence of a massive galaxy at the cluster cen- 
ter, fo llowing the method employed in Meneghett i et al.] 
( 2003 ) . The galaxy was mod eled with pseudo-isothermal 
mass distribution (see e.g. Donnarumma et al. 20111 
with a velocity dispersion of a = 300 km/s and a cut-off 
radius of 23 kpc. 

We use the deflection angle maps of the A383 SL model 
presented here, and use them to perform a ray-tracing 
simulation. We stress that such a simulation is com- 
pletely consistent with those performed for each simu- 
lated halo. A large number of artificial elliptical sources 
is used to populate the source plane at z s = 2, which 
are distributed on adaptive grids with increasing spatial 
resolution towards the caustics, in order to sample with 
greater accuracy the regions where sources are strongly 
magnified. By counting the sources that are lensed as 
giant arcs, we measure the lensing cross section, which is 
the area surrounding the caustics where sources must be 
located in order to produce images with length-to-width 
ratios larger than 7.5. The deflection angle maps also 
allow to measure the cluster median Einstein ring, de- 
fined as the median distance of the critical points from 
the cluster center. In the following discussion, we refer to 
the Einstein radius for sources at redshift z s — 2. Note 
also that the median Einstein radius is defined differ- 
ently than the simple effective radius of the critical area 
which is usually used and was implemented throughout 
this work in order to compare to other results. In this 
section only, we use the median Einstein radius, in order 
to be consistent with previous work based on these sim- 
ulations and since it us ually better correlates w ith the 
lensing cross-section (see Meneghetti et al. 2011). 

By doing this analysis, we find that A383 has a me- 
dian Einstein radius mc d = 19.0 ± 1" , and a lensing 
cross section a = 2.49±0.9x 10~ 3 h~ 2 Mpc 2 . Comparing 
these results with the distributions of 9 mc d and a of ha- 
los with similar mass in the MareNostrum Universe, 
we find that A383 is a remarkably strong lens, given its 
relatively small mass. Indeed, the majority of simulated 
clusters with 6 x 10 14 < M mr < 7 x 10 14 h~ l M® have 
much smaller critical curves and cross sections. For ex- 
ample, the mean Einstein radius and lensing cross section 
of such sample are 8.9" and 5.6 x 10~ 4 /i~ 2 Mpc 2 , respec- 
tively. The values measured for A383 exceed the maxi- 
mal values measured in the simulations, which are 18.14" 
and 2.26 x 10 -3 h~ 2 Mpc 2 , respectively. As an example, 
the distribution of medi an E instein radii fo r simulated 
cluste r s is s hown in Fig 14 As shown by Meneghetti 



et al. (2011), the lensing cross-section is tightly corre- 
lated to the median Einstein radius. We note also that 
each cluster in the MareNostrum Universe was pro- 
jected along three independent lines-of-sight in order to 
account for possible projection effects, i.e. for including 
cases where clusters are seen nearly along their major 
axis. The largest Einstein radii and cross sections are 
indeed produced by clusters whose major axis is almost 
perfectly aligned to the line-of-sight. Nevertheless, it is 
difficult to find a cluster that matches the mass and the 
strong lensing efficiency of A383 among the simulated 
halos. Extending the upper mass limit of the simulated 



CLASH Strong-Lensing Analysis of A383 



13 




E.med 



Fig. 14. — Distribution of the median Einstein radii of halos with 
6 X 10 14 h~ 1 M Q < M vir < 7 X 10 14 H^Mq extracted from the 
MareNostrum Universe snapshot at z ~ 0.19. The blue-shaded 
region shows the size of the median Einstein radius of A383 with its 
error-bar. See |5]for more details and the definition of the median 
Einstein radius. 

sample to 8 x 10 14 /i _1 M0, we find that only two cluster 
projections have Einstein radii and lensing cross sections 
larger than those measured for A383, which is still in the 
> 99% tail of the correspon ding distributions. 



Meneghetti et al. (2010b) showed that the concentra- 
tions estimated from the projected mass distribution of 
strong lensing clusters are on-average biased high, i.e. 
they are higher than the corresponding concentrations 
measured from the three-dimensional cluster mass distri- 



butions (see also Hennawi et al. 2007). Such bias depends 
on the cluster mass, redshitt, and lensing cross section. 
For a fixed mass and redshift, clusters with large lens- 
ing cross sections are typically affected by a more sever e 
concentration bias. As done by |Meneghetti et al. (2011), 
we use the MareNostrum Universe clusters to esti- 
mate a lower limit of the concentration bias for A383-like 
clusters. To do that, we select the numerically simulated 
halos with redshift and mass matching those of A383 and 
lensing cross section a > 10 -4 ft~ 2 Mpc 2 . This limit was 
set in order to have a statistically significant sample of 
simulated halos. The median ratio of C2D/C3D for these 
lenses is ~ 1.35. Thus, for objects with a lensing effi- 
ciency as high as in A383, the concentration measured 
from lensing is expected to be > 35% higher than their 
true 3D-concentr ation. This expectation agre es well with 
a recent work by Morandi fe Limousin | (p011|) estimatin 



the triaxial shape of A383. Morandi & Limousin ( 2011 
deduced by a joint analysis of X-ray and SL measure 
ments (which are commensurate with our analysis), a 
concentration of c V i r ~ 6.1, while we obtained ind eed a 
44% higher value in our 2D analysis, c V i r ~ 8.8 (see ^4.2 



6. SUMMARY 



In this work we have presented a new detailed lens- 
ing analysis of the galaxy cluster A383 in multi-band 
ACS /W FC3 images. Our well-established modeling 
method flBroadhurst et~aI1 |2005b| |Zitrin et aLl |2009b 



20101 |MerTen et al|2011HZitrin et al.|2011a|bTc) Tias iden- 
tified^ multiplv-lcnsed images and candidates, so thai 



ply-lensed images and candidates, so that 
in total 27 images of 9 different sources were incorpo- 
rated to fully constrain the fit. Though more lensed can- 
didates may generally be found in this lensing field with 



further careful effort, the resulting model is clearly fully 
constrained by these multiple systems. 

The accurate photometric redshifts of the newly found 
multiple-systems enabled by the extensive multi-band 
HST imaging allow for the most secure lensing-based de- 
termination of the inner mass profile of A383 to date, 
through the cosmological lensing-distance ratio, and im- 
ply a mass profile of dlogE/dlogr ~ —0.6 ± 0.1, simi- 
lar to other well-known relaxed clusters, and in excellent 
agreement with WL analysis from wide-field Subaru data 
(Umetsu et al., in preparation, see also Figure pi). In 
addition, we note that our mass profile is generally con- 
sistent with a recent joint l e nsing , X-ray ,and kinematic 
analysis by |Newman et aL (2011), out to at least twice 
the Einstein radius where our SL data apply. 

In Figure [3] we plotted the critical curves along with 
the multiple images found and used in this work. For a 
source at z s = 2.55, the effective Einstein radius — 
16.3 ± 2", or ~52 kpc at the redshift of the cluster. This 
critical curve encloses a projected mass of M — 2.4 ± 
0.2 x 10 13 MfT), in agreement with other publi shed results 
(e.g., |Smith et all2001[ |Newman etlD201 1 ) . 

We compared the properties ot A383 with clusters of 
similar mass drawn from the MareNostrum Universe 
numerical simulation (see Sj5j). We find that A383 is a 
remarkably strong lens, given its relatively small mass. 
The majority of simulated clusters 6 x 10 14 < M v i r < 
7 x 10 14 h^ 1 M Q have much smaller critical curves and 
lensing cross sections. The largest Einstein radii and 
cross sections are produced by clusters whose major axis 
is almost perfectly aligned to the line-of-sight. Even with 
this taken into account, it is difficult to find a cluster 
that matches the mass and the strong lensing efficiency 
of A383 among the simulated halos, so that A383 lies at 
the > 99% tail of the corresponding distributions (Figure 
14). Accordingly, for objects with a lensing efficiency 
as high as in A383, the concentration measured from 
lensing is expected to be > 35% higher than their true 
3D-concentration, in agreem ent with recent results (e.g., 
Morandi fc Limousin||201lj) . 

A383 is the first cluster observed and analyzed in the 
CLASH framework (see §1). As we have shown, despite 
the relatively small Einstein radius and correspondingly 
low number of multiply-lensed images, the remarkable 16 
filter imaging allowed us to immediately uncover several 
new multiple-systems. With a statistical sample of 25 
massive galaxy clusters being deeply imaged with HST, 
we should be able to test structure formation models with 
unprecedented precision. 
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